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I. INTRODUCTION 

Lattice calculations have made considerable progress 
toward a first-principles determination of the K — > tttt 
weak decay amplitudes [TJ [2] . The methodology is now 
in place, results for the 1 = 2 final state with a complete 
error budget are available [2], and complete results for 
the more challenging 1 = final states should become 
available in the next few years. At that stage we will 
finally learn whether and in what manner QCD can ex- 
plain the AI =1/2 rule and the observed CP-violation 
rate in K — ¥ tttt decays. 

Encouraged by this progress, it is natural to consider 
what information lattice calculations might eventually 
offer concerning the decays of heavier mesons. For ex- 
ample, the LHCb experiment recently reported evidence 
for CP-violation in (the difference of) D° — > 7r + 7r~ and 
D° — » K + K~ decays [3|. Although the rate is larger than 
naive expectations from the Standard Model (SM), there 
is, at present, sufficient uncertainty in the SM prediction 
for it to be consistent with the LHCb result (see, e.g. 
Refs. [4 - 9J ) . This raises the obvious question of whether 
a calculation using lattice methods is feasible. 

The aim of this paper is to take a first step in de- 
veloping the methodology for such a calculation. We 
show how, if one can ignore all but two-particle chan- 
nels, then a generalization of the work of Liischer, and of 
Lellouch and Liischer, would allow, in principle, a calcu- 
lation of the required matrix elements from lattice calcu- 
lations in a finite volume. In practice, however, channels 
with more than two particles are coupled by the strong 
interactions to 7T7r and KK, e.g. the four pion channel, 
and they cannot be ignored at center of mass (CM) ener- 
gies as high as the D°-meson mass (Mpo = 1865 McV). 
Thus our method would yield only semi-quantitative re- 
sults for the desired matrix elements. Nevertheless, it 
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is a necessary first step, and work is underway to extend 
the methodology to channels with multiple particles (see, 
e.g., Ref. US]). 

It is instructive to recall the three essential ingredients 
needed for the lattice calculation of K — > irir amplitudes. 
First, one needs to know the relation between the energies 
of two-pion states in a finite box and the infinite- volume 
scattering amplitude. This was worked out by Liischer 
in Refs. [TTHTl] (and generalized to a moving frame in 
Refs. [T5HT7] ). Second, one needs the relation between 
the matrix element that one can determine on the lat- 
tice, which connects a kaon to a finite-volume two pion 
state, and the infinite volume matrix clement which de- 
termines the decay rate. This was provided by Lellouch 
and Liischer in Ref. |18) (and generalized to a moving 
frame in Refs. [El Hi])- Finally, one must calculate the 
large number of Wick contractions that contribute, in- 
cluding several quark-disconnected contractions requir- 
ing special methods and high statistics. In this stage one 
also extrapolates to physical quark masses. This entire 
program has been carried out for the 1 = 2 final state [5] , 
and a successful pilot calculation has been done for the 
more challenging 1 = case pQ. 

The calculation of D° decays is considerably more chal- 
lenging. In particular, the first two of the three aforemen- 
tioned ingredients need to be generalized to account for 
the opening of many channels. If we focus on the 1 = 
final state, then strong-interaction rescattering connects 
two-pion final states to those with four, six, etc. pions, 
as well as KK and 7777 states^] As already noted, we 
consider here only the case in which several two-particle 
channels are open, which for the D° would mean keeping 



1 It is important to note that the fact that the D° has a very 
large number of decay channels [19] is not itself a concern, but 
rather that, having fixed the final-state quantum numbers, in our 
case to / = 0, there are still a large number of states. In a lattice 
calculation, one can separately consider the decays to states with 
differing strong-interaction quantum numbers. 
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the 7T7T, KK and rjrj channels while ignoring those with 
four or more pions. We make this approximation not be- 
cause we think that it is a good description of reality at 
the D° mass, but rather because it is a necessary first 
step towards the required formalism. 

Within this approximation, we provide here the gen- 
eralization of both the Liischer quantization condition 
and the Lellouch-Liischer (LL) formula. These general- 
izations are useful also in many other systems. For exam- 
ple, the quantization condition allows the determination 
of the parameters of the S-matrix in the 1 = channel 
above the two kaon threshold (and thus in the region of 
the / (980) resonance), because the coupling to four or 
more pions remains weak for such energies. The same 
should be true in the 1 = 1 case, where KK and ijtt 
are the dominant channels in the vicinity of the do (980). 
The multi-channel LL formula can be used to calculate 
K — > 7T7T amplitudes including isospin breaking (so that 
7T^~ and 7T° are not degenerate). Generalization to baryon 
decays are also possible, but this requires dealing with 
particles with spin, which we do not attempt here. 

There have been a large number of recent papers study- 
ing the generalization of the Liischer quantization condi- 
tion to multiple two-body channels [20 23 and assessing 
its utility. The work of Ref. [5U] uses non-relativistic 
quantum mechanics, while Ref. [21] is based on a non- 
relativistic effective field theory. References [55] and [53] 
are based on relativistic field theory, and give an explicit 
result [Eq. (3.5) of Ref. [22]] for the case of two s-wave 
channels in which the total momentum vanishes and in 
which the contributions from higher partial waves are as- 
sumed negligible. We also note that the multiple channel 
problem has been studied using an alternative approach 
based on the Bethe-Salpeter wavefunction 2.41. 

We provide here, as a step on the way to the gener- 
alized LL formula, a derivation of the multiple-channel 
quantization condition within quantum field theory. We 
include all allowed mixing between different partial 
waves. No assumptions about the form of the interac- 
tions are needed, aside from the proviso, common to all 
approaches, that the range of the interaction must be 
smaller than the box size. Also, our result holds for 
any value of the total momentum P of the two parti- 
cle system (i.e. it holds for a moving or a stationary 
frame), and for either degenerate or non-degenerate par- 
ticles in each channel. We follow closely the approach of 
Ref. [TH], which presented a generalization of Liischer's 
single-channel quantization condition to a moving frame. 
Indeed, we find that the most general form of the final 
result, given in Eq. (29 1, is identical in form to that of 



S-matrix which, due to unitarity and symmetry, is de- 
termined by three real parameters [see Eq.(48) below]. 



Ref. |16j (modulo some minor changes in notation). 

After deriving the general quantization condition in 
Sec. [Tl] we restrict our considerations to the simplified 
situation in which only s-wave scattering is included. We 
focus on the case with two channels (suggestively labeled 
7T7r and KK), although we also provide the generaliza- 
tion to more than two channels. In the infinite volume 
theory, the two-channel system is described by a 2 x 2 



We use a particular parametrization of S to rewrite our 
quantization condition in a convenient, pure real form 
[Eq. (49 1]. We explain how our result is equivalent to 
that of Ref. [22] in the case of a stationary frame. As 
is discussed in Refs. [51Tf53"] . three independent pieces of 
information are needed to determine the three indepen- 
dent S-matrix parameters at each center of mass energy, 
E* . References [22] and [23] discuss in some detail the 
prospects for using either twisted boundary conditions 
or uneven box sizes for this purpose. We restrict our- 
selves here to an alternative approach, also mentioned 
in Refs. [22] and [23], of using three different choices for 
the parameters {L, P}, where L is the box size. (We 
assume a cubic box and periodic boundary conditions.) 
The parameters {L, P} must be tuned such that there is 
a two-particle state in the spectrum having the desired 
value of E* . In this way one obtains three independent 
conditions, and can solve for the S-matrix parameters at 
the chosen value of E* . 

Turning now to the LL formula, we follow the same 
approach as used by Lellouch and Liischer in Ref. [18]. 
Specifically, we add a Z?-meson to our two-channel system 
and analyze the effect of an infinitesimal weak perturba- 
tion on the quantization condition. This yields a rela- 
tion between a finite-volume weak matrix element and 
a linear combination of the desired infinite-volume ma- 
trix elements. In Sec. |IV| we present a derivation of the 
relation which follows closely the original LL work. In 
the final result, Eq. (85), the coefficients relating finite 



and infinite volume matrix elements are given in terms 
of the S-matrix parameters and their derivatives, evalu- 
ated at the decay particle's mass. These can be calcu- 
lated using the multiple-channel quantization condition, 
as sketched above. It turns out that three different lattice 
matrix elements are needed to separately determine the 
two infinite- volume matrix elements. Note that this is the 
same as the number needed to determine the S'-matrix 
parameters. For more than two channels this correspon- 
dence no longer holds. 

Since the S'-matrix parameters and their derivatives 
are ultimately determined from the spectral energies, it 
should be possible to write a form for the generalized 
LL formula in terms of the spectral energies and their 
derivatives alone. We derive such a form in Sec. W\ The 



result, Eq. ( 122 ), is probably more useful in practice than 



Eq. (85). The second derivation also brings out an im- 
portant feature of the generalized LL formula. Finite- 
volume energy eigenstates in the coupled-channel theory 
can be written as linear combinations of infinite volume 
7T7T and KK states having (in our case) I = as well as 
the higher values of t allowed by the cubic symmetry of 
the box. The LL methodology is (as noted in the orig- 
inal paper) simply a trick to determine the coefficients 
of the relevant irir and KK states. This point has also 
been stressed recently by Ref. [55] in a different context. 
Our second derivation makes clear that, irrespective of 
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the details of the weak Hamiltonian, one always obtains 
the same linear combination of tttt and KK states, and 
that this feature holds for any number of channels. 

The remainder of this article is organized as follows. In 
the following section we give our derivation of the multi- 
ple channel quantization condition. In Sec. |III| we restrict 
to s-wave scattering and derive a useful form of the con- 
dition. The multiple channel generalization of the LL 
formula is then derived in Sec. |IV| and the alternative 
derivation is presented in Sec.|Vj We conclude in Sec. |VI| 
We include an appendix, in which we discuss the gener- 
alization of Watson's theorem to two channels. 

The generalization of Liischer's quantization formula 
to multiple channels for arbitrary P using field-theoretic 
methods (the work described in our Sees. [TT| and III I has 
also been considered by Briceno and Davoudi [26J. Our 
results are in complete agreement (although we use a 
different parametrization of the S-matrix). Their paper 
is being released simultaneously with the present article. 



II. MULTIPLE CHANNEL EXTENSION OF 
QUANTIZATION CONDITION 

In this section we derive an extension to multiple two- 
body channels of the Liischer quantization condition, 
which relates the infinite volume scattering amplitudes 
to finite volume energy levels. We assume throughout a 
cubic spatial volume with extent L and periodic bound- 
ary conditions. The (Minkowski) time direction is taken 
to be infinite. The total momentum 



ordering we need is 



P 



lirnp 
L 



(n P <= Z 3 ) 



(1) 



is fixed but arbitrary, i.e. the quantization condition we 
derive holds for a "moving frame" as well as a stationary 
frame. We first consider the case of only two open chan- 
nels, describing the extension to an arbitrary number of 
channels at the end of this section. 

We take each channel to contain two massive, spinless 
particles. The particles of channel one are labeled pi- 
ons and are taken to be identical with mass mj = M v . 
The particles of channel two, called kaons, are taken non- 
identical, though still degenerate, with mass m-i = Mk- 
What we have in mind is that the first channel corre- 
sponds to the I — tttt state, and the second to the 1 = 
KK state. Including both identical and non-identical 
pairs allows us to display the factors of 1/2 that appear 
in the former case. We consider degenerate particles to 
simplify the presentation, but describe the generalization 
to non-degenerate masses at the end of this section. 

For concreteness, and to match the physical ordering, 
we take the pion to be lighter than the kaon. For our 
results to hold, we must assume that the thresholds for 
three or more particles lie above the two kaon threshold. 
If we assume a G-parity like symmetry, so that only even 
numbers of pions can couple to a two-pion state, then the 



2M r < 2M K < E* < 4M ff , 
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where E* is the center of mass (CM) energy. The only 
possible scattering events are then 
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If E* drops below 2Mk , only the tttt channel is open and 
the problem reduces to that discussed by Liischer [TTHLl"] . 

The inequality 2Mr < 4M W does not, of course, hold 
for physical pions and kaons — the four and six pion 
thresholds occur below that for two kaons. Neverthe- 
less, the coupling to these higher multiplicity channels is 
weak at low energies, and our results should still hold ap- 
proximately as long as we are not too far above the two 
kaon threshold. Indeed, it may be that, in the 1 = case, 
the 7777 channel becomes important before that with four 
or more pions. If so, our formalism would still apply, 
generalized to three channels as described below. The 
approximation of ignoring channels with more than two 
particles will become increasingly poor as the energy in- 
creases, and will likely give only a rough guide by the 
D mass. A qualitative indication of this (ignoring dif- 
ferences in phase space) is that the /o(1500) has a 50% 
branching fraction to 4tt, while the branches to tttt, KK 
and 7777 are ~ 35%, 9% and 5%, respectively [T§] . 

The two channel quantization condition is obtained by 
a straightforward generalization of the single-channel ap- 
proach of Ref. [16] . To make this note somewhat indepen- 
dent of that reference, we reiterate some of the pertinent 
details. We begin by introducing a two body interpo- 
lating field <t(x) (not necessarily local) which couples to 
both channels. Following Ref. [TB] we then define 



cup) 



J(-P-x+Ex°), 



1a(xyt(0)|0) (4) 



where P = (E, P) is the total four momentum of the two 
particle system (in the frame where the finite volume 
condition is applied), and 



L ; x 



= d 



(5) 



is the spacetime integral over finite volume. The relation 
to the CM energy used above is 



E* 



E 2 -P 2 . 



(6) 



The poles of Cl give the energy spectrum of the finite 
volume theory, and thus the condition that Cl diverge is 
precisely the quantization condition we are after. 
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FIG. 1. (a) The initial series of ladder diagrams which builds 
up Cl [see Eq. Q]. The Bethe-Salpeter kernels iK are con- 
nected by fully dressed propagators. The dashed rectangle 
indicates finite volume momentum sum/integrals, (b) and 
(c) The series which build up the matrix element A and the 
scattering amplitude iM.. Note that these series contain only 
the momentum integrals appropriate to infinite volume, (d) 
The resulting series for the subtracted correlator [see Eqs. (171 
and (26 ll. Each dashed vertical line indicates an insertion of 



T, which carries the entire volume dependence (neglecting 
exponentially suppressed dependence). 



To proceed to a more useful form of the condition, 
we follow Ref. [IB] and write Cl in terms of the Bethe- 
Salpeter kernel, as illustrated in Fig.[TJa): 

+ f a j , q [z 2 A 2 ]. k JK M;q A^ A2 ]i m y m , q '+--- ■ 

J L\q,q' 



(7) 



The notation here is as follows. Indices j, k, I and m 
refer to the channel, and take the values 1 or 2. The two 
particle intermediate states are summed/integrated as is 



appropriate to finite volume 

-y 

L 3 



L,q 



2^ 



(8) 



The summand/integrand includes the product of two 
fully dressed propagators 

[z 2 A 2 ] i = S^zMMl)] [zi(P-q)Ai(P-q)] , (9) 



where 



A 3 -(«) 



q 2 — in 2 - + it 



(10) 
(11) 



Here <pi an d ^2 arc interpolating fields for pions and 
kaons, respectively, chosen such that Zj — 1 on shell. 
rji = 1/2 and 772 = 1 account for the symmetry factors of 
the diagrams. K is related to the Bethe-Salpeter kernel 



iK 



= tBS l] (q,P-q,-q , ,-P + q') 



(12) 



with BSij the sum of all amputated j — > i scattering dia- 
grams which are two-particle-irreducible in the s-channel 
(with particles of either type). Finally, a^ q and crt , de- 
scribe the coupling of the operators a and to the two- 
particle channel j. Their detailed form is not relevant; 
all we need to know is that they are regular functions of 

q- 

We emphasize two important features of Eq. Q. First, 
it does not rely on any choice of interactions between the 
pions and kaons, such as those predicted by chiral per- 
turbation theory. All the quantities that enter can be 
written in terms of non-perturbatively defined correla- 
tion functions. Second, the kernel iK and the propa- 
gator dressing function z have only exponentially sup- 
pressed dependence on the volume [12] . Thus, if L is 
large enough that such dependence is negligible (as we 
assume hereafter), we can take iK and z to have their 
infinite- volume forms. 

The dominant power-law volume dependence enters 
through the momentum sums in the two-particle loops. 
To extract this dependence, we use the identity derived in 
Ref. [16] , which relates these sums for a moving frame to 
infinite-volume momentum integrals plus a residue. Be- 
fore stating the identity we recall the relevant notation. 
For any four vector = [k , k) in the moving frame, 
= (k°*,k*) is the result of a boost to the CM frame. 
In particular, the total four-momentum (E, P) boosts to 
(E* , 0) in the CM frame. We also need the quantities 



q ; = ^(E*) 2 /4- m 2 



(13) 



which are the momenta of a pion (J = 1) or kaon [j = 2) 
in the CM frame. The identity then reads (no sum on i 
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here) : 



/(fc)i fc A i (fc)A i (P-%(A;) 



L:k 



+ j dn q ,dn q *, f*(f) Fu(q*,q*)g*(q*), (14) 



with 



/ -/ 

J oo;Aj J 



d 4 k 



(15) 



We introduce two functions f(k) and g(k) to correspond 
to the momentum dependence entering from the left and 
right of the loop integrals, as well as that from the dress- 
ing functions [see Fig. [TJa)]. The functions / and g must 
have ultraviolet behavior that renders the integral/sum 
convergent. In addition, the branch cuts they contain, 
corresponding to four or more intermediate particles, 
must be such that, after the k° contour integration, they 
introduce no singularities for real fc. This condition holds 
when < E* < AM^. The last line of (fl4|) depends 



on the values of the functions / and g when the two 
particles are on-shell, and thus only on the direction of 
the CM momentum, q* . Specifically, if is the moving 
frame momentum that boosts to the on-shell momentum 

(E*/2,qi*), then 



/*(?*) = /(%), 9t(q*)=9(qi) 



(16) 



Finally the quantity which depends on q* , L and the 
particle mass, contains the power-law finite-volume de- 
pendence of the loop sum/integralj^] Its form is given 
below in Eqs. (23)-(25). Note that it is diagonal in chan- 



nel space, i.e. it cannot change pions into kaons. It can, 
however, insert angular momentum, due to the breaking 
of rotation symmetry by the cubic box. 

The key point of the identity is that the difference be- 
tween finite and infinite volume integrals depends on on- 
shell values of the integrand, allowing the finite-volume 
dependence to be expressed in terms of physical quan- 
tities. Applying the identity to each loop integral in 
Fig. [lja), one then rearranges the series by grouping 
terms with the same number of insertions of T . The 
volume-independent term with no F insertions is of no 
interest, since it does not lead to poles. Thus we drop it 
and consider the difference 



c sub (p) = c l (p) - cup) ■ 



(17) 



In the remaining diagrams with F insertions, all terms 
to the left of the first F and to the right of the last are 
grouped and summed into new endcaps which we label 
Aj and Aj [see Fig. lib)]. These quantities equal certain 
matrix elements of the interpolating field a |16| 

Aj(k*) = (k*, -F; j; 0^(0)10),^, (18) 



4(**)=<0|«7(0)|fc*,-**; i; ^7^) |S , |=9 



(19) 



In contrast to [16| we include no wavefunction renormal- 
ization factors, because our single particle interpolating 
fields satisfy on-shell renormalization conditions. Having 
summed up the ends the next step is to do the same for 
the series which appears between adjacent F insertions 
[Fig. [ljc)]. As indicated in the figure, this series gener- 
ates the infinite volume scattering amplitude iAiij. We 
thus deduce an alternative series for C su i, built from A, 
A' and iMs, all connected by .Fs [Fig. [TJd)]. 

We stress that the analysis just performed is a straight- 
forward generalization of the single channel analysis of 
Ref. [H] . All that has changed is that F and M. are now 
2x2 matrices in channel space, and A and A' vectors. 

To proceed, we decompose A, A' , M and F in spherical 
harmonics, defining coefficients via 



(20) 
(21) 

(22) 



Aj(k') = V47rA i;Am y^ m (r) 
A'jCk*) = V^A'^Yl^k*) 



Mij{k*,k*) 



4rtMij ; e 1<mi; t 2 , m2 Y ei!mi (k*)Y? (k *) 



J~ij(k , k ) — 



1 % , ^ (23) 



where a sum over all Ps and m's is implicit. The factors 
of 47r are present so that we match the conventions of 
Ref. [16]. They imply, for example, that for a purely 
s-wave amplitude, M. is the same in the two bases (for 
the 4-7T cancels with the two spherical harmonics). The 
kinematical factor F is given in Ref. [16 (aside from the 
above-noted factor of rji) and takes the forrrrl 



Fij;£ 1 ^rni;£ 2 ,m 2 — &ij Pi;£i ,mi \l 2 ,m 2 



8nE* 1- ^ 2 miTO 2 



i 

2E* 



WE4? c Lte 2 ) / dnY; um Yi m Y i2 



(24) 



2 The result {wj is equivalent to Eqs. (41-42) of Ref. [16], although 
we have done some further manipulations to the last line of | |14[ l 
to bring it into a matrix form. Also, we have included a factor 
of rji in J 7 , rather than keeping it explicitly as in Ref. |16| . 



3 An additional difference from Ref. [16| is the appearance of Re q* 
rather than q* . This is discussed in the next section. 
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Here the volume-dependence enters through the sum^] 



uil e 



*(q* 2 -k> 2 ) 



L 3 ' ujk q* 2 — k* 2 
k 



d 3 k* 



k*W4nY e , m (k*) 



'>/(> 



-k' 2 ) 



(2tt) 3 q* 2 -k* 2 



(25) 



with Wfc = y k 2 + m 2 being the energy of a particle with 

momentum k, and u>l the energy after boosting to the 
CM frame. The properties of these sums are discussed in 
Ref. [S]. 

We are now in a position to write down the final result. 
The series indicated in Fig. [ljd) gives 



C auh (P) = -J2 A ' F HMF] n A, 



n=0 



-A'- 



F- 



iM 



-.A. 



(26) 



(27) 



Here all indices are left implicit and may be restored in 
the obvious way. For example, 



A'FMFA — 4(i,mi- f, i];(l,mi;(2,m ! 



Mi 



jk]&2 >77l2 ;^3 ,1713 -^fcZ;^3 ,m3 ;^4,m4 -^-1,^4 ,7714 



A, 



(28) 



As Coo has no poles in the region of E* that we consider 
(below &M n ), the poles in Cx must match the poles in 
C su b- The desired quantization condition is then just that 
the matrix between A' and A have a divergent eigenvalue. 
This may be written as 



det (F- 



iM) = 0. 



(29) 



where we recall that the matrices now act in the product 
of the two-dimensional channel space and the infinite- 
dimensional angular-momentum space. More precisely, 
F is diagonal in channel space but has off-diagonal ele- 
ments between different angular momentum sectors (as 
allowed by the symmetries of the cubic box and the mo- 
mentum P), while M is diagonal in angular momentum 
but off-diagonal in channel space. 

Equation (29 1 is the main result of this section. It 
has exactly the same form as that for the single channel 
given in Ref. |16j (aside from the change of notation in 
which symmetry factors are contained in F rather than 
kept explicit). The generalization to more than two two- 
particle channels is now immediate. As long as E* is kept 
below the four particle threshold of the lightest particle 
the arguments above go through in the same manner. 



One need only extend the values of the channel indices, 
taking care to include the appropriate symmetry factor 
r/j for each channel. The final result then has exactly the 
form of Eq. p9l. 



To make the formal expression ( 29 ) useful in practice 
that there 



one assumes that there is some £ ma x 
partial wave amplitudes are negligible 



M 



= 0. 



above which the 



(30) 



One can then show that, although F couples £ < £ max 
to £ > £ max , the projection contained in M is sufficient 
to collapse the required determinant to that in the £ < 
l max subspace. The argument for this result is given for 
one channel in Ref. [H] and generalizes trivially to the 
multiple channel case. Thus one finds that Eq. (291 still 
holds, but with M. and F now understood to be finite 
dimensional matrices both in channel space and in the 
partial wave basis, with £ running up to £ max - 

To conclude this section we comment briefly on two 
generalizations of the result. We first consider the case 
when not just a single a but rather a set of operators {cr a } 
is of interest. This is likely to be the case in practice since 
multiple operators may be needed to find combinations 
with good overlaps with the finite-volume eigenstates. If 
there are n such operators, then Cl generalizes to an 
n x n matrix: 



Cf(P) = 



a i(-P-x+Ex°) , 



\o*{x)J\0)\Q) ■ (31) 



L:x 



The generalization of Eq. (27) is effected by replacing A' 
with an n x 2 matrix 



(^'i A' 2 ] 



I A[ a=1 
V i 



'a=i\ 



(32) 



and A with a 2 x n matrix 



A 2 



A\= l AI 



b=2 



A b 2 = 2 



(33) 



The key point, however, is that the matrix between A' 
and A is unchanged, so that the quantization condition 
(291 is unaffected. This is as expected, since the opera- 



tors used to couple to states cannot affect the eigenstates 
themselves. 

The second generalization is to the case of non- 
degenerate particles. The expressions given above re- 
main valid as long as one makes three changes. First, the 
symmetry factors rji become unity for all non-degenerate 
channels. Second, q* in Eqs. (24 1 is replaced by the so- 
lution of 



We are slightly abusing the notation here for the sake of clarity. 
c^ m depends not only on q* 2 but also on m;, but we keep the 
latter dependence implicit. The dependence is made explicit at 
the end of this section. 



E* 



M 2 



M 2 

lvl ib i 



(34) 



which is the CM momentum when the channel contains 
particles of masses Mi a and M^. Third, when evaluating 
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c im using Eq. ( 25 ), one should use one of the masses M ia where 



or Mib when determining ojk, and k* . One can show 
that both choices lead to the same result. 

The third change emphasizes that the kinematic func- 
tions cf m depend not only on q* but also on the particle 
masses. This can be made explicit by rewriting them in 
terms of a generalization of the zeta-function introduced 
in Ref. [H]. The result is 



cL(g* 2 ) 



4?r / 2tt 

r e Y em (f) 



fV2 



(r 2 — x 2 ) s 



Zf m [l;(q*L/2ir) 2 }, (35) 
(36) 



where 7 = E/E* , n runs over integer vectors, and r is 
obtained from n by 7*11 = 7 _1 [ri|| — cnp] and r± = n±, 
where parallel and perpendicular are relative to P, and 
2c = (1 + (M 2 a - M 2 b )/E* 2 ). The sum in Z lm can be 
regulated by taking s > (3 + 1)/2 and then analytically 
continuing to s = 1. This shows that mass dependence 
enters through the differencer] M 2 a — M 2 b . One can derive 
( 35 1 by generalizing the method used for the degenerate 
case in Ref. ITBl. 



III. MULTIPLE-CHANNEL QUANTIZATION 
CONDITION FOR S-WAVE SCATTERING 



For the remainder of this article we focus on the sim- 
plest case, £ m ax = 0, in which only s-wave scattering is 
significant. In this section we determine the explicit form 
for the finite-volume quantization condition when there 
are two channels. We also present compact forms for 
the condition when an arbitrary number of two particle 
channels are open. 

With only s-wave scattering, the two channel quanti- 
zation condition takes the form 



■iMi 



[iM s 12 ][iM s 21 ]=0, (37) 



5 The apparent lack of symmetry under the interchange M; a -o- 
Mih can be understood as follows. One can show that Zf — > 
under this interchange (so that for degenerate masses 
the zeta-functions for odd I vanish 1151 ). This sign flip for odd 
I must hold also for the cf , and it does because the inter- 
change of masses leads to k* — > — k* at the pole. The sign flip 
is canceled in the expression for T , Eq. | |23[ ), since the product 
^fijmi m2 (k *) a ' so changes sign. This is because, when 

I is odd, the integral over d£l in the definition of _F, Eq. {24}, 
enforces that t\ +I2 is odd. The overall effect is that the quanti- 
zation condition is symmetric under mass interchange, as it must 
be. 



Reg* 



c P (q. 



8irE* 

2\ 



I 

2E* 



1 00 



(9i 



(38) 
(39) 



and the superscript on F and M is a reminder that only 
1 = contributes. 

To simplify Eq. (37), and in particular to re-express 



it as an equation between real quantities, it is useful to 
recall first the single-channel analysis. This has the ad- 
ditional benefit of showing how the two-channel result 
collapses to the known single-channel result in the ap- 
propriate kinematic regime, namely 



2M W < E* < 2M K . 



(40) 



In this regime q 2 becomes imaginary, and th e se cond 
channel contributes negligibly because c p [Eq. (25)] be- 
come s ex ponentially volume-suppressed and Re q* in F 2 
[Eq. (24)] vanishes]^] Sending F 2 — > wc find that the 
quantization condition becomes 



[M' u ]- l = m 



SirE* 



2E 



(41) 



Note that the pion momentum q\ is real for the energy 
region considered. 

Naively one might think that Eq. ( 41 ) gives two con- 



ditions, the separate vanishing of the real and imaginary 
parts. This is not the case, however, because the vanish- 
ing of the imaginary part is a volume-independent condi- 
tion which is guaranteed to hold by the unitarity of the 
5-matrix. This can be seen by expressing M. in terms of 
the real phase shift S(q*), 



8ttE* 



e 2i<5( 9 *) _ ^ 
2i 



V1Q1 
8ttE* 



[cot <J fa?) -i] 



. (42) 

Here e 2ld is the one dimensional unitary S'-matrix in the 



partial wave basis. Given Eq. (42 1, it is manifest that the 
imaginary part of Eq. (41) holds automatically. The real 



part of ( 41 ) then gives the moving frame generalization of 



the Luscher result in the familiar partial wave form [T31 

main] 



tan<5(<?i) = — tan(/> p (g*) . 



where 



tan 4> p (q*) 



4?r 



r p i * 
c (q 



(43) 



(44) 



6 The appearance of Re q* rather than q* in Fi can be understood 
by reviewing the derivation of F in Ref. I16| . The term enters as 
the difference between principal part and ie prescriptions. When 
q* is imaginary there is no pole and different ways of regulating 
give the same result. 
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We now return to the CM energies for which both 
channels are open, 2Mk < E* < &M n , and generalize 
Eq. (43 1. The first step is to recall the relationship be- 
tween the scattering amplitude and the S-matrix. Uni- 
tarity implies that 



M 



iM s ^P 2 M s 



(45) 



where P 2 is a diagonal matrix containing the phase-space 
factors, whose square root is 



P = 



(46) 



We note that, when expressed in terms of q* , the form 
of P is still valid if the two particles in the channel are 
non-degenerate. We also note that the form (45) holds 



for an arbitrary number of two-particle, s-wave channels, 
with P generalized in the obvious way. 
The solution to the unitarity relation is 



iM s = P- 1 (S s - 1) P~ 



(47) 



where 5 s is a dimensionless, unitary 2x2 matrix. To 
proceed, we need to parametrize S s . First we note that 
S s can be taken to be symmetric. This is because of the 
T-invariance of the strong interactions, together with the 
fact that angular momentum eigenstates have definite T- 
parity (in our case, positive). Thus in the 2x2 case, 
S s is determined by three real parameters. We use the 
"eigenphase convention" of Blatt and Biedenharn [30 , 



S s = 



e 2i5 a o 
e 2iS ' 3 



c e s ( 
-s £ c. 



(48) 



where the notation s^ = sin x and c^ = cos x will be 
used throughout. The three real parameters 8 a , 8g, and 
e generalize the single 8 which appears in the one channel 
case. The parameter e quantifies the mixing between the 
mass eigenstates of channels one and two (the pions and 
kaons) and the S'-matrix eigenstates. The phases S a and 
Sg are, for arbitrary e, associated with both channels. Of 
course, in the limit e — > they reduce, respectively, to 
the phase shifts of pion and kaon elastic scattering. 

Substituting the form of S s into Eq. ( j47j ) and then 
placing this in Eq. ( 37 ) and simplifying, we deduc^] 



[tan<5 a + tan(/) P (q , 1 )] [tan Sg + tan <t> p (q^)] 

+ sin 2 e [tan5 Q — tan<5^] [tan</) P (g^) — tan<^> p (<72)] 

= . (49) 



7 We emphasize that the physical content of Eq. j49| , namely that 
there is a relation between scattering amplitudes and energy lev- 
els, does not depend on the para metr ization chose n for the matrix 
S a . This is clear either from Eq. ( |37| ) or from Eq. | |58[ ) below. An 
adva nta ge of our choice of parametrization is that it shows that 
Eq. |37j only implies one real condition (rather than two), an 
observation which must hold for any parameterization. We also 
note that the freedom to independently change the phases of ■kit 
and K K states, which leads to S s — > U^S S U, with U a diagonal 
unitary matrix, does not change the quantization condition, as 
can be seen most easily from Eq. 1581 below. 



This is the main result of this section. One can use it in 
one of two ways: to predict the spectrum given knowledge 
of the scattering amplitude from experiment or a model, 
or to determine the S'-matrix parameters from a lattice 
calculation of the spectrum. In the former case, we note 
that all quantities appearing in (49), i.e. 5 a , Sg, e, q* 
and 4> p , are functions of E* , One can thus search, at 
given spatial extent L and total momentum P, for values 
of E* which satisfy Eq. ( 49 ) . If the condition holds for a 
particular Ef, then 



E k (L;n P ) = 



P< 



(50) 



is in the spectrum of the finite-volume moving-frame 
Hamiltonian. Here we choose to write E k as a function 
of Hp rather than P, since, in practice, it is the former 
quantity which is held fixed as one varies L. 



The second use of (49) is the most relevant for the 
discussion in subsequent sections. For a given choice of 
E* , one finds, through a lattice calculation, three pairs 
{L,np} for which there is a spectral line Pfc such that 
El [defined in Eq. (50)] is equal to E* . One could use 



a fixed Hp and consider multiple spectral lines (the sim- 
plest choice conceptually), or use three different choices 
of np (probably more practical since one would not need 
to determine so many excited levels). In cither case, one 
ends up with three versions of Eq. ( 49 ) , all containing the 



desired quantities S a (E*), Sg(E*) and e(E*), but having 
different values of the 4> p (q*). Solving these equations 
one determines, rather than the angles themselves, the 
quantities tan<5 a , tan Sg, and sin 2 e at CM energy E* . 
For our discussion we therefore restrict 



S a ,B G [-7r/2,7r/2] and e G [0,tt/2] 



(51) 



Having determined the restricted phases over a range of 
energies, one can afterward relax the constraints in order 
to build continuous curves as a function of energy. We 
direct the reader to Refs. [53] for discussion of other 
methods for extracting the three scattering parameters. 

We emphasize that Eq. (49) has a very intuitive form. 
If S a = Sg or mi = m 2 or t = then the second line 
vanishes and the result reduces to two copies of the one 
channel quantization condition [Eq. (43)]. To see that 
this makes sense, note that for identical phase shifts, the 
e matrix commutes through the phase matrix and we 
recover two uncoupled channels. Similarly if the masses 
are degenerate then the eigenstates of the S-matrix will 
also be mass eigenstates leading to the decoupled form. 
Finally, the decoupling for e = is an obvious property 
of the parametrization [Eq. (48)]. 



An alternative solution to the unitarity relation (45) 
can be given in terms of the _ftT-matrix used in Ref. |22) . 
Specifically, (|45| is satisfied if 



{M'r 1 = M -iP 2 /2, (52) 
with M any real symmetric 2x2 matrix. If we set 
1 



M 



8ttP* 



Iv^ll^llv^ll 



(53) 
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[where double bars denote a diagonal matrix, so that 
|^|| = diag(y / ?yi, y/rfe)], and further set P = 0, then 
one can show that the two-channel quantization condi- 
tion given above can be manipulated into the form given 
in Eq. (3.5) of Ref. (55] in terms of the real, symmetric 
matrix K. 



We now generalize Eq. (49) to N s-wave channels. As 



noted above, the form of the unitarity relation (45 1 holds 
for any N, and the same is true for its solution (47). In 



the latter, the N channel S'-matrix can be parametrized 
as0 



S = Rr 



2iS I 



with R an SO(N) matrix, and 



„2iS\ 



: diag[ 



2iS r 



R. 



2iS p 



(54) 



(55) 



Together with Eqs. (54) and (55) one needs the N x N 



generalization of F, which has been discussed in the pre- 
vious section. From these definitions one can straightfor- 
wardly work out the quantization condition for N cou- 
pled channels. 

We conclude this section by describing two additional 
ways of writing the quantization condition, both of which 
make the higher channel generalization especially clear. 
Observe first that, for any number of open channels, 



F- 1 = P- 1 111 



P- 



Combining this with (47), it follows that 



F~ 



iM = P~ 



S 



P- 



(56) 



(57) 



Since P _1 has no singularities in the kinematic regime 
we consider, the quantization condition can be rewritten 

as 



det S 



-2i<t 



0. 



(58) 



This form shows that the symmetry factors cancel from 



the quantization condition in general. Although Eq. (58) 



looks like it will lead to one complex and thus two real 
conditions, it turns out that it leads only to a single real 
condition. This follows from the identity 



1 + iU 



S- 



-2i<t> I 



2l 



iT^l - it s \\R 
R 1 II \\R + \\td> II 



(59) 



where t x = tanx. It gives a manifestly real rewriting of 
the quantization condition, 



det R' 1 iff \\R+ iJ = 



(60) 



The remainder of this paper is limited to the s-wave, so we drop 
the superscript s hereafter. 



This form leads directly to the result (49) in the two 



channel case, and collapses to the single-channel result 
(43) for any channel that decouples from the rest. If 



any of the channels contain non-degenerate particles, this 
enters only through the values of the kinematic functions 
t«A, as discussed in the previous section. 



IV. MULTIPLE CHANNEL EXTENSION OF 
THE LELLOUCH-LUSCHER FORMULA 

Having found the two channel quantization condition, 
we are now in a position to work out the two channel gen- 
eralization of the LL formula which relates weak matrix 
elements in finite and infinite volume. The derivation fol- 
lows the original work by Lellouch and Liischer, Ref. [18] , 
which was extended to a moving frame by Refs. [H)I I17j . 

We begin by introducing a third channel which is de- 
coupled from the original two. This contains a single 
particle, which we call a Z3-meson, whose mass satisfies 



M D > 2M 7r ,2M K 



(61) 



We next introduce a weak perturbation to the Hamilto- 
nian density 



H(x) H(x) + \H w (x) , 



(62) 



where A is a parameter which can be varied freely and 
can, in particular, be taken arbitrarily small. The per- 
turbation Jiw is defined to couple channels one and two 
(pions and kaons) to the third (D-meson) and nothing 
more. It is convenient to choose it to be invariant under 
time reversal (T) symmetry. The generalization to per- 
turbations which are not T invariant is described at the 
end of the section. 

Consider now the finite volume spectrum, first in the 
absence of the perturbation. The spectrum of two- 
particle states with P — 2i:np/L is determined by 
Eq. (49). It is L-dependent and L can therefore be tuned 



to make one of the levels, call it kn, degenerate with the 
energy of a single (moving) D meson 



E kT 



M%+P 2 



(63) 



With L fixed in this way, we now turn on the weak in- 
teraction. At leading order in degenerate perturbation 
theory this changes the energies to 



E {i) =eW ±\V\M W \ 



(64) 



where V — L 3 , E^ ' = Ek D , and the finite- volume matrix 
element is 



M w = L (k D \H w (0)\D) L . 



(65) 



The subscripts L on the states indicate that they are nor- 
malized to unity, unlike the relativistically normalized in- 
finite volume states. Superscripts W are used throughout 
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iAM lk = 



FIG. 2. The diagram giving rise to the amplitude perturba- 
tion AM [See Eq. |[69|], 



this section to indicate that the quantity includes both 
the leading order and the 0(A) correction, while super- 
scripts indicate the unperturbed quantity. The effect 
of the perturbation may also be written in terms of the 
CM energy as 



E* {1) = M D ± XAE* 
E^V\M W \ 



AE* = 



Md 



(66) 
(67) 



Of course, in addition to affecting the finite volume 
energy spectrum, the weak perturbation also changes the 
infinite volume scattering amplitudes. The leading order 
effect is generated by the diagram of Fig. [2] which gives 



M {1) =M {0) tAATW 



with 



AMj, k = 



(j\H w (0)\D)(D\Hw(0)\k) 



2E(°)V\M, 



w 



(68) 



(69) 



The indices j and k run over the two two-particle chan- 
nels. This perturbation may equivalently be represented 
through shifts in S a , 5p and e 



=6^(E*)±\A6 a (E*) 



5^(E*) = S^(E*)±XASp(E*) 



r(0) 



£ W(ET) = e w {E*)±XAe(E* 



(70) 
(71) 
(72) 



The explicit forms of the perturbed phases are given in 
Eqs. pl|-(pl) below. 



The calculation now proceeds as follows. When the 
quantities 



8$\E*), 6f(E*), ande(°)(£*) 



(73) 



are placed in the quantization condition [Eq. (49)], the 
condition is satisfied by construction at E*^ 3 ' — Mjj. 
Alternatively if one places 



(74) 



into the same condition then it must be satisfied when 
evaluated at the perturbed CM energy E*^\ but only 
to linear order in A. The constant order term in the A 
expansion is just the unperturbed condition, and so it is 



the vanishing of the 0(A) term that is of interest. The 
condition that this term vanish gives the relation between 
finite and infinite volume weak matrix elements that we 
are after. 

The only detail left to discuss, before substituting into 
the quantization condition and expanding in A, is the ex- 
plicit forms of the amplitude corrections to S a , Sp and 
e. Before these are found it is useful to determine the 
constraints on the infinite volume matrix elements which 
arise from Watson's theorem. As shown in App. [Aj time 
reversal invariance and unitarity constrain the matrix el- 
ements to be such that the following two quantities are 
real: 



V2 



-iSp 



q\ri\A D ^-K c e + q 2 r] 2 A D ^ KK s e , (75) 

- V^lAo-Wr S e +y/ q\j] 2 A D ^ KK C e . 



Here 



.4 



D— >7T7T 



\H W (0)\D), 



(76) 



(77) 



and similarly for the KK case, normalized so that the 
decay rates to each channel are 



8ttM2 



1 



2M T 



D- 



(78) 



This relation holds also if the particles in a channel are 
non-degenerate (requiring rj — 1). All energy-dependent 



parameters in (75) and (76), i.e. S a , Sp, e and q* , are to 
be evaluated at E* = M n . 



The results ( 75 ) and ( 76 1 hold when the phases of the 



states are chosen so that the S'-matrix is symmetric (as 
is possible given T invariance). This does not determine 
the signs of the two matrix elements, and these signs are 
unphysical. More precisely, the relative sign ambiguity 
is the same as the ambiguity in the sign of e, so once we 
have fixed the latter to be positive, the relative sign is 
physical. The overall sign remains unphysical, and can 
be chosen, for example to set v\ > 0, 



Inverting the relations ( 75 ) and ( 76 1 yields 



Ad^it = ~jL= he* c £ -v 2 e l5 ? s e ] (79) 



Ad^kk = r— [vie 1 "- s e +v 2 e 1 ^ c e 



(80) 



Inserting these in AA4 , Eq. ( 69 1 , and using the relation 



between Ai and S, Eq. ( 47 1 , and the parametrization of 



S, Eq. (48 ), we find that perturbations to 8 a , dp and e are 
real. This is a consistency check on the calculation (or 
an alternative derivation of the Watson's theorem con- 
straint). Specifically, we find 



A5 a = 
ASp = 



-Mvt 



V\V 2 



Cq cp(t a —tp) 



(81) 
(82) 

(83) 
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where t Q = tani5 Q , etc., and 



16irE(°)M D V\Mi 



(84) 



We now have all the ingredients to substitute into the 
quantization condition and determine the LL generaliza- 
tion. We emphasize that, when the expansion in A is 
performed, 5 a , Sp and e each contribute both from the 
amplitude corrections of Eqs. (|8ll)-( 83 ) and from the shift 



in the energy, ( 66 1 . The other contributions arise from 
the energy dependence of <\>i = (j) P (q*). Substituting and 
simplifying, we find the main result of this section 



where 

Cl2 

Cm 2 



C\v\ + C 2 v 2 +C\2V\V2 = C M 2 \M W \ 2 

7T tl+t2+2t,3+(t 2 -tl)(l-2s^ 

16 cl 

tt t 1 +t 2 +2t a +(ti-t 2 )(l-2sg) 

16 c 2 

4 C a Cp 

ir 2 M D V 2 {E^) 2 [ti 



(85) 

(86) 
(87) 



* (t 2 +t (3 +(t Q -t^ S 2) 
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+ 4 (ti+t Q +(t^-t Q )sf) 
+ ^(t 2 +t^+(ti-t 2 )^) 



+ 



4t^ 
M D 
4 s?' 



M 



D 



ti+t a +(t 2 -ti)s?J 

(tl - t 2 )(t a — tg) , 



and where we use the additional notation 



tan i 



(89) 



(90) 



All quantities are evaluated at the D mass, and we have 
dropped the superscript The primes on 4>i indicate 
derivatives with respect to q* while those on <5 a , Sp and 
e indicate derivatives with respect to E* . In each case, 
these are the natural variables on which the quantities 
depend. We have checked that this formula reduces to 
(two copies of) the single-channel LL result if e — > 0. 



We now describe how the result (85) can be used in 
practice. A lattice calculation yields the finite-volume 
matrix element |Afy^|, and the aim is to determine the 
infinite- volume matrix elements Au^-k-k an( i 4^^. 



Using the generalized quantization condition ( 49 1 for 



three different spectral lines (all chosen to have E* — 
Md) one can determine 5 a , 5p and e as described in the 
previous section. Repeating the procedure at a slightly 
different energy allows a numerical determination of the 



required derivatives. One now evaluates \M\y \ at the de- 
generate point on one of the spectral lines. The knowl- 
edge of the S'-matrix parameters and their derivatives, 
together with the value of L, allows one to calculate the 
values of the four C's [Eqs. ( 86 )-(89 )]. Combined with the 
value of I M\y \ , one then finds from Eq. ( 85 ) a quadratic 
constraint on v\ and i> 2 . Repeating the procedure for 
a second spectral line gives an independent constraint, 
which allows for the determination of v\ and u 2 up to 
a two-fold ambiguity corresponding to the unknown rel- 
ative sign. Finally, repeating for a third spectral line 
resolves the sign ambiguity. With v\ and w 2 determined 
in this way, one can obtain the infinite-volume matrix 
elements using Eqs. (79) and (80). Although this pro- 



cedure is rather elaborate, we note that (for the case of 
two channels) three spectral lines are needed both for the 
determination of the parameters of the S'-matrix and of 
the LL factors. 



We conclude this section by commenting that Eq. ( 85 ) 
factors as 



where 



sgn(Ci)(ci«i + c 2 v 2 ) 2 = C M A M w\ 



ci = vjCij c 2 = sgn[CiCi 



(91) 



(92) 



The only new information encoded in Eqs. ( |9l[ ) and (92) 
relative to Eq. ( 85 ) is that 



4CiC 2 = C 



12 ! 



(93) 



which can be shown to hold by applying Eq. ( 49 ) to 



Eqs. (86)-(88). Although the factorized form (91) is sim 



pier, it does not reduce the number of values of L that 
are needed because there remains a sign ambiguity (from 
the square root) at each L. What it does make clear, 
however, is that the generalized LL condition will fail 
when the signs of C\ and Cm 2 are opposite. Presumably 
this cannot happen for physical values of the phase shifts. 
We stress that this issue also arises in the original one- 
channel set-up, where the LL formula only makes sense 
if 



d(S + f_ 
dq* 



> 0. 



(94) 



We return to these sign constraints in the next section. 
The form (91) also allows one to write the LL condi- 



tion as a factored form in terms of the desired matrix 
elements, 



\c w A 



c K A 



D-+KK\ 



\Cm 2 \\M w \ 



(95) 



where and ck are complex, and can be determined 
from the above results. As this equality holds for any 
T-invariant form of weak perturbation and for any de- 
cay particle, it must imply a relation between finite and 
infinite volume states 



,{ko \ oc 0^(^71,0^1 + ck{KK, out I 



(96) 
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Here the ellipsis indicates the hit and K K states of higher 
angular momentum which are needed to satisfy the peri- 
odic boundary conditions. Indeed, as noted in the orig- 
inal derivation of Ref. [TS], the use of the D-meson is 
simply a trick to obtain the normalization factors]^] It 
follows that Eq. (95) must also hold for perturbations 
which are not T- invariant. 



The appearance of the linear combination in Eq. ( 96 ) 



can be better understood from an alternative derivation 
of the LL formula, to which we now turn. 



we have 



X^ x =0. 



(101) 



Now, while holding L fixed, we change the energy by 
±XAE = ±\V\M W \ and change M to M {Q) T XAM, 
and require that the quantization condition still hold. 
Thus we have, to linear order in A, 



det(# (0) + XAX) = 0, 



(102) 



where 



V. ALTERNATIVE DERIVATION OF 
LELLOUCH-LUSCHER FORMULA 

In this section we present a different derivation of the 
two channel LL relation which has the following advan- 
tages: (a) it does not require determining the shifts AS a , 
Ad/3 and Ae, but rather works directly with the change in 
M. ; (b) it gives one directly a condition with the factored 
form, proportional to the left hand side of Eq. (95); (c) it 



allows one to rewrite the LL condition in a simpler form 
in which the only inputs required are the derivatives of 
the energies with respect to L along the three spectral 
lines. This form is likely to be more practical. 

We work directly with the condition det(F _1 + = 
0, and keep results for general number of channels, N, as 
far as possible. We begin by defining 



X = F~ L 

y = s- 



iM 

,-2#| 



and recall from Eq. (57) that 



x = p^yp- 1 



(97) 
(98) 



(99) 



The quantization condition det X = is equivalent to 
X having an eigenvector with vanishing eigenvalue. We 
label this eigenvector e x . Note also that the symmetry 
of X implies (~e? x ) Tr = '^e x is a left eigenvector, also 
with zero eigenvalue. 

Now we can formulate the LL condition in a relatively 
compact form. As above, let be the scattering am- 

plitude at CM energy E* = M D . Similarly, let be 
the finite-volume factor at this CM energy and for one 
of the values of box size L for which the quantization 
condition holds. Then for 



*(°) = (F(o))-i +iM (o) 



(100) 



AX = ±AE 



dX 
dE 



T- iAM . 



(103) 



It follows that there must be a new eigenvector of the 
form 



l? x 



XAlt x 



(104) 



which is annihilated by the perturbed matrix. From the 
0(A) term in 



\<F X 



\At x ] 



X [0) 



AAA" 

we deduce 

^ x AXlt x =0. 
Using the explicit form of AX this becomes 
dX 



XAlt x ] =0, 
(105) 

(106) 



AE ^ x 



dE 



^ x 



= Y x iAM ~S 



(107) 



where the derivative is evaluated at E* = Mp. 

We can slightly simplify this result by expressing the 
left hand side in terms of y rather than X, and thus 
removing factors of The point is that, when the 

derivative acts on the P -1 factors in X, the contribu- 
tion to the left hand side vanishes, since one can still act 
(cither to the left or the right) on the zero-eigenvector. 
Thus we can rewrite the condition in terms of the zero 
eigenvector of y, which is 



P-^ X 



The new form is 



AE^ 9y 



dE 



<F X 



[iAM]d 



(108) 



(109) 



9 In the one-channel case, an alternative line of argument has been 
developed for obtaining the LL relation, based on matching the 
density of two-particle states in finite and infinite volumes |31| . 
In the present case, we do not see how to use this approach to 
determine the relative normalization, cx/cvr, of the two compo- 
nents in the finite volume state. Thus we think that this ap- 
proach could provide only a consistency check. 



We now focus on the 2x2 case. To proceed, we need 
the explicit form for ~e? Y , which is given, up to an overall 
normalization factor, by 



^ Y = 



1 



(110) 
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where z is the real quantity 



z = t f 



sin 



+ 4>x) 



sin(5 / 3 + 02 ) 
It is clear from the form of Eq 



(|109D 

Y 



(111) 



and the relation 



( 108 1 that the normalization of ~e ry is irrelevant and so 
we have chosen a relatively simple unnormalized form. 

We evaluate the right hand side of Eq. ( 109 ) using the 
form of AM., 



Eq. 



( 69 1 . It is immediately apparent that 
the result factorizes, given that AM is an outer product. 
This will hold for all N. In the N = 2 case we have 



< e x [iAM}~e 
where 

Moo = 



-2i0i 



M 



8nM D E(°)V\M w \ 



(112) 



■KK 



(113) 



Here we have used the assumed T-invariance of Hw- We 
have pulled out the phase e -2 ^ 1 so that is real. 

Its reality is not obvious, but can be established using 
the results derived from Watson's theorem and given in 
App. [A] In particular, an algebraic exercise shows that 



sin 



-vi 



i '2 - 



sin(6 Q + 2 ) sin(^+(j6 2 ) 



(114) 



and we recall that the quantities V\ and V2 [defined above 
in Eqs. (75)-(76l] are real. 



The result (113) makes clear that, for any choice of 
Hw, one ends up with the matrix element to a given 
(complex) linear combination of (irir\ and (KK\ states, 
since all the factors {<f>\ 1 <p2 and z) are determined by 
and L. Indeed, what the LL method has allowed 
us to do is determine the coefficients of the s-wave (mr\ 
and (KK\ components within the finite-volume state. As 
mentioned above, this decomposition has nothing to do 
with Hw, and thus we can use the result for any Hw, 
including one involving T-violation. By comparing the 
result (113) to the general decomposition of the finite- 
volume state, Eq. (96), we can read off the ratio of the 
coefficients, 



ck 

(-it 



2-0l) 



(115) 



It is interesting that the relative phase between ck and 
c w is determined by the kinematic phas es <pj . Given the 
form of AM, and the fact that, in Eq. (|112[), it is sand- 



wiched between e x and ~et x , it follows that the zero 
eigenvector itself gives the relative size of the irir and 
KK contributions: 



CK 



(116) 



This illustrates in a direct way that the linear combina- 
tion which appears is completely independent of the form 
of Hw, since the eigenvector of X knows nothing about 
this perturbation. 

Having discussed the right hand side of Eq. ( 109[ ) in 
some detail we now turn to the left. Specifically, we show 
that it is possible to write the left-hand side in terms of 
the derivative of the spectral energy with respect to L. 
To motivate this form, we first recall that the LL re- 
sult of the previous section depends on S Q , 8» and e and 
their derivatives, evaluated at E* = Mo- As described 
in Sec. |IV| the three ^-matrix parameters may be de- 



termined, using Eq. (49), by finding three different pairs 
{L, ftp} for which there is a spectral line E k (L;fip) sat- 
isfying Ef.{L\np) — Md [see Eq. (50)]. Furthermore, by 



slightly changing the three L values, one can determine 
S a , 5/3 and e at slightly different energies and thus deduce 
the derivatives at Mp,. 

The point of reiterating these steps is to note that, 
since the lattice simulation actually gives the energy spec- 
trum as a function of L, it would be preferable if the LL 
result could be rewritten directly in terms of the prop- 
erties of the spectrum. In this way the extra step of 
separately working out the phase shifts and their deriva- 
tives would be avoided. This turns out to be possible, as 
we now show. 

We use the quantization condition in the form det y = 
0. To stay on a spectral line E^{L]np) as we vary E 
away from the moving frame D-meson energy E^°\ we 
need to vary L in such a way that this condition remains 
fulfilled. We note that, while F depends on both E* and 
L, S depends only on E* . Thus we use E* and L as 
independent variables. Then the condition to stay on a 
spectral line becomes 



= t y 



which leads to 







. oy 




dE* 


L 


^ ~d~L 





dEl 
dL 



line 



^Ydy-^Y 

c dL c 



1?r 



(117) 



(118) 



Here, in the left-hand side, the subscript "line" indicates 
that the derivative is along a spectral line with fixed ftp. 
The key features of Eq. ( 118 ) are that the denominator 



on the right-hand side is, up to a simple overall factor, 
equal to the quantity appearing on the left hand side of 
the Eq. (109), while the numerator is a kinematic factor. 



Specifically, using 



oy 

dE* 



^ydy^y 

E dE • 



and 



dEk 
dL 



E* k dE* k 



line 



E k dL 



line 



EL 



(119) 



(120) 
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(which follows since E 2 = (E* ) 2 + (P L) 2 / L 2 and because 
PL is fixed along the spectral line), we find 



dE 



-2ie 



-2i0i 



d<t>i , ,2 d4> 2 
dL~^__dL 



dE k I 
dL I line 



_p2_ 

EL 



Combining this with (109) and (112) we conclude 
M 2 



dL 



16irM D E(°)V 2 \M w \ 2 
We thus have found an alternative form of the LL relation 



T Z dL 
dE k I , 
dL I line ^ £(°)L 



(121) 



(122) 



which is simpler than Eq. (85 1, and also likely to be more 
practical. 



The single channel version of (122) is instructive. It 



can be written, using Eq. ( 78 1 , in terms of the decay 
rate: 



D- 



2E^V 2 \M W \ 

Mn 



dL 



dEk_ I , P 2 
dL I line B(.°)L 



(123) 



This form holds both for identical and non-identical par- 
ticles, with the symmetry factor being contained in T. It 
also sheds light on the sign constraints discussed in the 
previous section. The right-hand side must be positive. 
Based on numerical studies, we find that d(j)/dL is al- 
ways positive, implying that the denominator, which is 
proportional to dE^/dL, must be negative. 



The same holds for the two-channel result, Eq. (122) 



In order for the right-hand side to be positive, the de- 
nominator must be negative. Since we could do the LL 
analysis on almost any spectral line, this appears to im- 
ply that dEt/dL < in general. The only exception 
is for a state with E* k below the two particle threshold. 
Such a state occurs, for example, as the lowest energy 
state for P = if there is an attractive interaction. For 
such a state one has dE* k jdL — dEk/dL > 0, i.e. of the 
"wrong" sign. But in this case the LL analysis does not 
apply, because the particle lies below threshold in infinite 
volume. 



VI. CONCLUSIONS 

We have presented two new results in this article. 
First, a field-theoretic derivation of the generalization of 
Liischer's quantization formula to the case of multiple 
strongly coupled two-particle channels (where the par- 
ticles are spinless). Second, the generalization to mul- 
tiple channels of the Lellouch-Liischer formula relating 
finite- volume and infinite- volume matrix elements. We 
also have explained in some detail how, in the case of 
two channels, one can use these results to determine the 
infinite-volume decay amplitudes of a particle which is 
coupled by a weak interaction to the two-body channels. 

As already noted in the introduction, this is but a step 
on the way toward our "dream" application, namely the 
calculation of D° — > itit and D° —> KK amplitudes. To 



achieve that goal, one will also need to include the chan- 
nels with four or more pions. These are significant once 
one approaches the energy Mq. Work in this direction 
is underway. 

An example where our formalism should be useful 
with minimal approximation is the determination of the 
isospin breaking in K — > irn decays. Given the mass 
splitting between charged and neutral pions, there are re- 
ally two two-body channels to consider, and in this case 
the coupling to the four pion channel is very small and 
can reasonably be neglected. 
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Appendix A: Two-channel Watson's theorem 

In this appendix we work out the consequences of Wat- 
son's theorem for the phases of the matrix elements of in- 
terest, (mT\n w (0)\D} and (KK\U w (0)\D). We assume 
at first that T-Lw is T invariant, and describe the gener- 
alization to non-invariant Hamiltonians at the end. We 
closely follow the textbook presentation given in Ref. [35] . 

We consider the 3x3 S-matrix with the three states 
being the hypothetical D meson (at rest) and the s-wave 
7T7T and KK states. We assume that we are in the kine- 
matic regime described in the main text, so that the 3x3 
S'-matrix is unitary. Although we introduce a weak cou- 
pling between the D and the two particle states, so that 
the D is a resonance, its width is of second-order in the 
weak interaction and thus can be ignored at the linear 
order to which we work. Thus it is valid to treat it as an 
asymptotic state. 

Watson's theorem follows by breaking the S'-matrix 
into a strong part S<°) and a weak part S w . The strong 
part is T invariant, and, since we use states which have 
definite (positive) T-parity, can be taken to be symmet- 
ric. This fixes the phases of the 7T7t and KK states, 
though not their overall signs. Extending the dimen- 



sionless, strong-coupling S-matrix of Eq. (47) to include 
the D gives 



S (0) = 



1 

s s 



(Al) 



where 1 is the lxl identity and S s is the 2x2 s-wave 
S-matrix given in ( 47 1 . The weak part only contains 



couplings between the D and the two-particle states, and 
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qW 

7T7r,D 
qW 

D KK,D 



aW qW 

D D : TTTT °D,KK 



(A2) 



we find 



w 

jk^>k,D 



iR-ikS, 



2iSi 



(A8) 



The assumed T invariance implies that it, too, is sym- 
metric. The non-zero elements of S w are proportional 
to the desired matrix elements 



s™ D = c p jj {j\[-m w (o)}\D) 



(A3) 



where j = 1,2 runs over the itit and KK channels, P 
is the square root of the phase space factor defined in 



Eq. (46), and c is a known real constant whose value will 



not be needed. 

Unitarity of the complete ^-matrix implies that the 
terms linear in the weak interaction satisfy 



iS w = S W(iS w ) i SW. 



This implies that 



k,D) i 



(A4) 



(A5) 



where in the last step we have used the symmetry of S n 
Using the explicit form for the two-channel /S-matrbP" 



5 s = R~ 







R, 



(A6) 



It follows that the phase of iRjkSYo is e 4,5j . This is 
the desired generalization of Watson's theorem to two 
channels. Thus the quantities 



-V4ttE* iRjkS^D 



(A9) 



are real. Using ( A3 ) we can rewrite the Vj as in Eqs. ( 75 ) 
and (1761). 



If the weak interaction is not T invariant, then Sj^ D will 
contain some number of T-violating phases. Since we are 
working to linear order in the weak interaction, wc can 
break up T-Lw into parts each with a single T-violating 
phase and treat each separately. Each such part has an 
overall phase e l< ^ r , and the symmetry of the S-matrix is 
replaced by 



Sr>.k( 



qW 

D k.D 



(A10) 



with 



R = 



-S f C, 



(A7) 



However, if we first pull out the overall phase by hand, 
then the symmetry of AS is restored, and Watson's the- 
orem applies to the residue. 
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